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ABSTRACT 

Gravitational instability has been invoked as a possible mechanism of giant planet formation in pro- 
toplanetary disks. Here we critically revise its viability by noting that for the direct production of 
giant planets it is not enough for protoplanetary disks to be gravitationally unstable. They must 
also be able to cool efficiently (on a timescale comparable to the local disk orbital period) to allow 
the formation of bound clumps by fragmentation. Combination of dynamical and thermal constraints 
puts very stringent lower limits on the surface density and temperature of disks capable of fragment- 
ing into self-gravitating objects: for the gravitational instability to form giant planets at 10 AU gas 
temperature at this location must exceed 10 3 K for a minimum disk mass of 0.7 Mq and minimum 
disk luminosity of 40 L Q . Although these requirements relax in more distant parts of the disk, masses 
of bound objects formed as a result of instability even at 100 AU are too large (~ 10 Mj) to explain 
characteristics of known extrasolar giant planets. Such protoplanetary disks (and planets formed in 
them) have very unusual observational properties and this severely constrains the possibility of giant 
planet formation by direct gravitational instability. 

Subject headings: planets and satellites: formation — solar system: formation — planetary systems: 
protoplanetary disks 



1. INTRODUCTION. 

Recent discoveries of Jupiter-like planets around solar- 
type stars have rejuvenated the interest in the issue of 
the origin of giant planets. Core instability model (Pcrri 
& Cameron 1974; Harris 1978; Mizuno 1980) in which 
Jupiter-like planets acquire their massive gaseous atmo- 
spheres by unstable gas accretion onto the preexisting 
massive solid cores has been one of the most fruitful ideas 
in this field. For rather long time this avenue of planet 
formation did not seem compatible with the short ob- 
served lifetimes (10 6 — 10 7 yr) of protoplanetary disks be- 
cause of the long time needed for the core accumulation. 
However, recent work by Rafikov (2003) and Goldreich 
et al. (2004) has found core formation time to be actu- 
ally quite short, essentially removing the timescale issue 
from the core instability scenario. Nevertheless, gravi- 
tational instability (hereafter GI) in the protoplanetary 
disk (Cameron 1978; Boss 1998) has been put forward as 
an alternative mechanism of giant planet formation. In 
this model massive gaseous disk becomes gravitationally 
unstable and rapidly fragments into a number of self- 
gravitating bound structures, which further collapse to 
become giant planets. A number of recent hydrodynam- 
ical simulations (Boss 1998; Mayer et al. 2002) employing 
isothermal equation of state (EOS) have confirmed this 
general picture and provided ample support for the idea 
of planet formation by GI. The goal of this study is to 
constrain this avenue of planet formation by putting spe- 
cial emphasis on the conditions necessary for the actual 
fragmentation of protoplanetary disk into self-gravitating 
objects. 

2. DYNAMICAL AND THERMAL CONSTRAINTS. 

Gravitational instability in a Keplerian disk can oper- 
ate only when the gas sound speed c s , surface density £ 
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and local angular frequency £1 satisfy the condition 



D.c 



ttGE 



<Qo- 



(1) 



(Safronov 1960; Toomre 1964). Here Q is the so-called 
"Toomre Q"; sound speed is defined as c s = (kT/fi) 1 / 2 , 
where T is the disk midplane temperature, fj, is the gas 
molecular weight, and k is a Boltzmann constant. An- 
alytical arguments and results of numerical simulations 
suggest that Qq w 1 . Equivalent way of formulating the 
condition for operation of GI is to require 
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where p is a midplane gas density and o,au = a /(l AU) 
is the distance from the central star scaled by 1 AU. 

Dynamical constraint Q is a necessary condition for 
GI to set in. However, even when QJ is fulfilled, gi- 
ant planet formation becomes possible only if the disk 
can actually fragment into bound self-gravitating objects. 
Recent studies (Gammie 2001; Rice et al. 2003) have 
demonstrated that fragmentation is possible only pro- 
vided that the cooling time of the disk t coo i satisfies 
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where £ is a parameter of the order of unity; numerical 
simulations (Gammie 2001; Rice et al. 2003) suggest that 

Cooling time of the disk can be estimated as the ratio 
of its thermal energy to the escaping radiative flux: 
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where a is the Stephan-Boltzmann constant, T is the 
midplane disk temperature, r rs kE/2 is the disk optical 
depth (k is the opacity), and 7 is the adiabatic index of 



gas. Function /(r) describes the efficiency of disk cool- 
ing and its specific form in (@J corresponds to the case 
when cooling is radiative. When the disk is optically 
thick, t ^> 1, radiation has to leak out through the large 
optical depth of the disk material; this lowers the effec- 
tive temperature at the disk photosphere by a factor of 
T 1 ' 4 compared to the midplane temperature T and makes 
tcooi very long [first term in the definition of /(r)]. In 
the optically thin case, r < 1, according to the Kirch- 
hoff's law (Rybicki & Lightman 1979) disk emissivity is 
low and t coo i is again very large [second term in the def- 
inition of /(r)]. Thus, disk cools most effectively when 
r»L What is most important for our further discussion 
is that the value of /(r) is above unity for any r and any 
cooling mechanism because effective temperature of the 
disk cannot exceed its midplane temperature. 

Dynamical constraint (JJJ limits the value of the sound 
speed from above. At the same time, expressing the tem- 
perature in through c s we find that J3J sets a lower 
limit on the gas sound speed. Combination of these con- 
ditions leads to the following constraint on c s necessary 
for the giant planet formation by GI: 
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where ( = £(7 — 1) rs 1. Only when (JSJ) is fulfilled 
can the disk be gravitationally unstable and cooling be 
fast enough for fragmentation to allow the formation 
of bound gaseous clumps, which later collapse to be- 
come giant planets. Similar argument has been advanced 
by Levin (2003) in application to self-gravitating disks 
around AGNs. 

At a specific location in the protoplanetary disk the 
condition Q can be satisfied only provided that the gas 
surface density obeys 
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Here fl = fx/mn is the mean molecular weight relative 
to the atomic hydrogen mass m# and M* is the mass of 
the central star (M@ is the Solar mass). For molecular 
gas of solar composition p, an 2.3 and T, min rj 3.4 x 10 5 g 
cm" 2 at 1 AU [for /(r)/(Q^ 4 C) = 1], while if hydrogen 
is atomic /2 rs 1.2 and X m m w 5.6 x 10 5 g cm -2 . In 
both cases the disk surface density at 1 AU turns out to 
exceed that of the minimum mass Solar nebula (Hayashi 
1981) by more than 10 2 , implying that a very massive 
protoplanetary disk is necessary to sustain giant planet 
formation by GI. 

According to (JjjJ, whenever (JJJ) is fulfilled, the sound 
speed in the disk is also bounded from below: 
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Because of that midplane temperature has to satisfy 

T>T mm = T mf [f(T)fl b where (9) 
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In the form given by (|6*|l- (|10J) constraints on the disk 
surface density and temperature still depend on the be- 
havior of opacity. However, bearing in mind that /(t) > 
1 for any r, one immediately sees that values of £j n / and 
Ti n f represent absolute opacity-independent lower lim- 
its on the disk surface density and temperature. Thus, 
we conclude that for the gravitationally unstable disk to 
be able to fragment (which is necessary for giant plan- 
ets formation) disk has to satisfy at least £ > £j n / and 
T > Ti n f. These lower limits are extremely robust since 
they do not depend on a specific mechanism of energy 
losses from the disk — by radiation or by convection. In 
practice, taking into account even the most basic proper- 
ties of the energy transfer within the disk one can formu- 
late much more stringent constraints, as we demonstrate 
below. 

3. APPLICATION TO PROTOPLANETARY DISKS. 

Protoplanetary disk luminosity Ld and mass Md are 
the characteristics which can be directly observed (or 
constrained using observations). According to JJjJ, mass 
of a protoplanetary disk which is capable of planet for- 
mation in the range of semimajor axes from o,„ to a out 
has to satisfy 
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This and all subsequent numerical estimates assume 
Qo = 1 and ( = 1. Note a rather weak dependence 
of Md on the disk dimensions. For the inner disk cutoff 
at Oi n = 0.1 AU and outer disk cutoff at a out = 100 AU 
one finds that Md > 3 M Q (for jj, = 1). This estimate as- 
sumes maximum cooling efficiency throughout the whole 
disk, i.e /(r) = 1, but even such lowest possible limit on 
Md well exceeds the mass of the central star. 

This constraint on the disk mass can be somewhat re- 
laxed if disk is producing planets only locally. Indeed, the 
typical scale length of the fastest-growing perturbation 
in the marginally gravitationally unstable disk (Q ~ 1) is 
A m 2irh, where h = c s /Cl is the vertical disk scaleheight. 
Thus, to form planets by GI it is, in principle, sufficient 
that only an annulus of the disk with the radial width of 
order A has surface density and temperature exceeding 
E m j n and T m i„. Using JSJ we can estimate 
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Because of the high midplane temperature (necessary to 
ensure efficient cooling) disk is not very thin geometri- 
cally; as a result, GI is likely to have global character. 

1 From infimum — the greatest lower bound of a set. 



Mass of such a "minimum planet-forming annulus" cen- 
tered on a and having width A is constrained by 



Af a >O.8Af a^ /5 [/(r)] 2/5 r 8/5 



mA 2/5 , * 

Mi) - ,13) 

Clearly, this lower limit on M a depends only very weakly 
on the semimajor axis of annulus. It also amounts to 
a sizable fraction of M± [even for molecular gas with 
(X = 2.3 and /(r) = 1 one finds M a > 0.2 M Q at 1 AU]. 
At the same time, typical masses of protoplanetary disks 
inferred from mm and infrared observations vary between 
1CT 3 M Q and 0.1 M© (Kitamura et al. 2002), but these 
masses of gas are typically extended over ~ 100 AU in 
disk radius (and not within just a narrow annulus). Be- 
sides, in the real nebula such annulus can easily be either 
very optically thick or very optically thin, which trans- 
lates into large value of /(r), additionally increasing the 
lower limit on M a (see below). 

Energy flux emitted from both sides by a unit surface 
area of the disk is given by 
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where we made use of (flH|l . Since dLd/dS oc a~ 24 / 5 
most of the energy is radiated by the innermost part of 
the gravitationally unstable disk. Total luminosity of a 
disk with an inner cutoff at a is 
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This estimate of luminosity holds even if only an annu- 
lus of gas is considered instead of the full disk, which 
is a direct consequence of the very steep dependence of 
dLd/dS on the distance from the central star. 

We now consider what these limits imply for the disk 
properties at different locations in the protoplanetary 
nebula. 

3.1. Limits on the disk properties at 1 AU. 

First we consider the possibility of giant planet forma- 
tion by GI in the region of terrestrial planets, at a = 1 
AU. Constraints obtained in previous sections imply that 
planetary genesis at this location requires rather ex- 
treme disk properties: temperature has to exceed at least 
Tint = 5.2 x 10 3 K (this estimate is based using fit, = 1.2 in 
|10j since gas cannot be molecular at such temperature) , 
surface density must be above Sj„/ = 5.7 x 10 5 g cm~ 2 , 
and disk luminosity has to exceed 10 4 L Q (all assuming 
C = 1 and Qo = 1). Depending on whether an extended 
disk or just a minimum size annulus is considered, min- 
imum mass varies from Md ~ 1.5 M© to M a rs 0.6 M Q . 
Apparently, even these least radical requirements are in 
complete disagreement with the observed properties of 
protoplanetary disks. 

The argument can be significantly sharpened by notic- 
ing that the disk with such high surface density has to be 
optically very thick, meaning that £ m j„ and T min pro- 
vide much more stringent limits on the disk properties 
than T,i n f and Ti n f. Using opacity dependence n{p, T) on 



the temperature T and gas density p (given by 0) from 
Bell & Lin (1994) we can substitute r = k(p, T)E m j n /2 
into © and © and solve the resulting system for £,„„ 
and T m i n . Performing this procedure one finds that at 1 
AU disk has to be extremely hot so that opacity is due to 
electron scattering. Disk temperature has to exceed 10 6 
K, but this is impossible, because such disk would not 
be bound to the central star and its radiation pressure 
would far dominate over the gas pressure. This firmly 
rules out the possibility of giant planet formation by GI 
within several AU from the central star. 

3.2. Limits on the disk properties at 10 AU. 

At 10 AU, in the region of giant planets, temperature 
will probably be low enough for the gas to be molecular. 
Using p, — 2.3 we find Ej„/ = 2.7 x 10 3 g cm -3 and 
Ti n f — 220 K. Disk luminosity has to exceed only 3.4 
L©, and the full disk and minimum mass annulus have 
to contain at least M4 = 0.4 M Q (for a out = 100 AU) and 
M a = 0.13 Mq correspondingly. These limits are more 
reasonable than at 1 AU, although mass is still high and 
disk is too hot compared to the observed systems, which 
typically have M d < 0.1 M Q and T < 10 2 K at 10 AU 
(e.g. Kitamura et al. 2002). 

However, using again opacities from Bell & Lin (1994) 
and p as given by J5J) we find that real disks (i.e. cooling 
not at the maximum efficiency) should be considerably 
more extreme. There are several possible solutions for 
Tmin and E ? „i n at 10 AU corresponding to different opac- 
ity regimes. The least extreme one is a "cold" solution 
with r a 60, T mm « 1100 K, S mm a 6 x 10 3 g cm" 2 , 
Ld > 40 L Q , M a > 0.7 M© which corresponds to molec- 
ular gas at the temperature of grain evaporation (other 
solutions have T m i n -J 7 x 10 3 K). This again yields a 
minimum disk which is too massive and hot to satisfy 
current observational constraints. Thus, it is extremely 
unlikely that GI can allow disk fragmentation and sub- 
sequent giant planet formation even at 10 AU. 

3.3. Limits on the disk properties at 100 AU. 

We also look at the possibility of giant planet formation 
by GI in the distant regions of protoplanetary nebula, at 
a = 100 AU. For molecular disk at this location we find 
E m/ = 20 g cm' 3 , T inf = 14 K, L d > 5 x 10' 3 L Q , and 
M a > 0.08 M Q . These properties change only a little 
when inefficiency of disk cooling is properly accounted 
for: using k rs 0.1(T/10 K) cm 2 g _1 in agreement with 
observations of protoplanetary nebulae (Beckwith et al. 
1990; Kitamura et al. 2002) we find that to be able to 
form giant planets by GI disk has possess at least r pa 2 
(marginally optically thick), T m i n ~ 20 K, S m i„ « 25 
g cm" 2 , L d > 10~ 2 L©, and M a > 0.1 M Q at 100 AU. 
Although M a is very near the upper end of the observed 
distribution of protoplanetary disk masses, these param- 
eters seem to be acceptable from the observational point 
of view. There are however additional reasons to doubt 
the possibility of planet formation by GI even at 100 AU, 
which we discuss next. 

3.4. Fragment masses. 

Another important observational constraint on the 
planet formation by GI comes from comparing the ob- 
served masses of extrasolar giant planets with the typi- 
cal masses of fragments into which disk breaks up when 



tcooi < £fi -1 . As we mentioned earlier, the lengthscale 
of the most unstable mode is A « 2irh, which results in 
the minimum fragment mass of 
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At 100 AU molecular disk meeting the requirements 
outlined in H3 3.3I would fragment into self-gravitating 
clumps with the mass of roughly 5 Mj[/(t)] 3 / 5 , where 
Mj is the Jupiter's mass; at smaller semimajor axes Mf 
would be higher. Even for /(t) = 1 this mass is larger 
than masses of most extrasolar giant planets detected to 
date (Marcy et al. 2003) . More realistic cooling efficiency 
corresponding to the optical depth of r sw 2 at 100 AU 
leads to Mf rs 9 Mj, landing minimum fragment mass 
not too far from the brown dwarf regime. 

It is possible that such clumps would be able to further 
fragment into smaller objects but we view this outcome 
as rather unlikely. Indeed, as the fragment contracts its 
optical depth increases and cooling time in this nonlinear 
regime becomes larger than the dynamical time of col- 
lapsing fragment, which disfavors subsequent fragmen- 
tation (Goodman & Tan 2004). It remains to be seen 
how the centrifugal support of the collapsing clump can 
change this conclusion. 

4. DISCUSSION. 

Novel analytical constraints presented in SJ21 an d El 
when confronted with observations of protoplanetary 
disks, severely undermine the possibility of giant planet 
formation by GI. In particular, we have demonstrated 
that disks capable of producing giant planets by GI at 
a distance of several AU from the central star cannot 
exist simply on dynamical grounds — to cool efficiently 
they must be too hot to be bound to the central object. 
This essentially rules out a possibility of in situ forma- 
tion of close-in extrasolar giant planets (" hot Jupiters" ) 
by GI. Rafikov (2004) have previously presented argu- 
ments against in situ formation of "hot Jupiters" via the 
core instability. Then the most natural way to explain 
the existence of close-in giant planets is to accept that 
they have formed elsewhere under more favorable condi- 
tions and then migrated to their current locations. 

Planet formation by GI is also extremely unlikely 
within several tens of AU, as results of ^ 13 3.21 demon- 
strate: disks with required properties must be so hot 



(T > 10 3 K), luminous (several tens of £©), and massive 
(~ Mq) that they would clearly stand out in a sample of 
observed protoplanetary nebulae. Beyond about 100 AU 
minimum disk properties allowing planets to form by GI 
become roughly acceptable from the observational point 
of view, although disk masses still reside at the very up- 
per end of the observed distribution of protoplanetary 
nebulae masses. However, it is still very unlikely that 
the extrasolar giant planets that we see now could have 
been produced by GI even at several hundred AU from 
the central star. The problem is not only in a poten- 
tial difficulty of migrating such planets from beyond 100 
AU all the way in to several AU, but also in producing 
planets with the right mass. As our estimate l|16f) of the 
minimum mass of unstable disk fragments demonstrates, 
bound objects produced by GI even at 100 AU are too 
massive (~ 10 Mj) to explain the observed mass distri- 
bution of extrasolar planets. 

Our study emphasizes the importance of the proper 
treatment of disk thermodynamics (especially its cool- 
ing) for studying the possibility of Jupiter-like planet 
formation. By now virtually all simulations which were 
able to demonstrate disk fragmentation and collapse of 
resulting dense objects in gravitationally unstable disks 
used isothermal EOS (e.g. Mayer et al. 2002). How- 
ever, use of this EOS is equivalent to setting the disk 
cooling time to zero which artificially relaxes the require- 
ments for the planet production process and is mislead- 
ing. Not surprisingly, more realistic calculations follow- 
ing thermodynamics in greater detail typically do not 
exhibit fragmentation of gravitationally unstable disks 
which are not capable of cooling efficiently (Pickett et 
al. 1998; Gammie 2001; Rice et al. 2003). Thus, simu- 
lations employing isothermal EOS should not be trusted 
too much when planet formation in real protoplanetary 
disks is concerned. 

Future infrared and mm observation will show whether 
protoplanetary disks with extreme properties satisfying 
the constraints necessary for giant planet formation by 
GI really exist. 
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